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Abstract 

It is shown in the framework of the operator product expansion and 
the renormalization group method that the twist-3 part of flavour nonsin- 
glet spin structure function Q2(x, Q 2 ) obeys a simple Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equation in the large Nc limit even 
in the case of massive quarks (Nc is the number of colours). There are 
four different types of twist-3 operators which contribute to includ- 
ing quark-mass-dependent operators and the ones proportional to the 
equation of motion. They are not all independent but are constrained 
by one relation. A new choice of the independent operator bases leads 
to a simple form of the evolution equation for g 2 at large Nc- 
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In the experiments of the polarized deep inelastic leptoproduction we can obtain 
the information on spin structures of nucleon, which are described by the two func- 
tions gi(x,Q 2 ) and g 2 (x,Q 2 ). The QCD effects on gi and g 2 have been extensively 
studied since earlier papers @-@|. Increasingly accurate measurements of g± 
have been performed at SLAC, CERN and DESY ||, while the g 2 measurements 
still have limited statistical precision ||. 

In the language of the operator-product-expansion (OPE) the twist-2 operators 
contribute to g\ in the leading order of 1/Q 2 - As for the structure function g 2 , on 
the other hand, both twist-2 and twist-3 operators participate in the leading order. 
Moreover, the number of participating twist-3 operators grows with spin (moment 
of g 2 ). Due to increase of the number of operators and the mixing among these 
operators the analysis of the twist-3 part of g 2 turns out to be rather complicated [[?[]- 
[ |H[| . In other words, the Q 2 evolution equation for the moments of the twist-3 part 
of g 2 cannot be written in a simple form but in a sum of terms, the number of which 
increases with spin. 

For the case of the twist-3 flavour nonsinglet g 2 it has been observed by Ali, Braun 
and Hiller (ABH) that in the large Nq limit g 2 obeys a simple Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equation | 16| . In their formalism of working di- 
rectly with the nonlocal operator contributing to the twist-3 part of g 2 , they showed 
that local operators involving gluons effectively decouple from evolution equation 
for large Nq- In fact their analysis has been made with massless quarks. 

In this paper I reanalyze the Q 2 evolution of the flavour nonsinglet twist-3 part 
of g 2 in the framework of the standard OPE and the renormalization group (RG) 
with massive quarks. Actually the OPE analysis of g 2 has been performed already 
and the anomalous dimensions of the relevant twist-3 operators have been calcu- 
lated PI P| |jn| [ fT/H [|T8| . However, to the best of my knowledge, the large N c limit of 



g 2 has not been thoroughly studied so far in OPE and RG. There are four different 
types of twist-3 operators which contribute to g 2 , including quark- mass-dependent 
operators and the ones proportional to the equation of motion. They are not all 
independent but are constrained by one relation. It was pointed out recently by 
Kodaira, Uematsu and Yasui [17] that any choice of the independent operator bases 
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leads to a unique prediction for the moments. Taking a new basis of the independent 
operators, I will show that the Q 2 evolution of the twist-3 part of g 2 obeys a simple 
DGLAP equation in the Nc — ► oo limit and thus the ABH result on g 2 is reproduced 
even with massive quarks. 

The spin structure function g 2 receives contributions from both twist-2 and twist- 
3 operators. However, the twist-2 part of g 2 can be extracted once g± is mea- 



!, 2 "- 2 U.Q- 2 ): - !h U.Q 2 )+ I 9 -^^dy. (1) 



sured [f[9| : 



Thus the difference 

g 2 (x,Q 2 )=92(x,Q 2 )-g t 2 w - 2 (x,Q 2 ) (2) 

contains the twist-3 contributions only. 

The twist-3 operators which enter the OPE for the flavour nonsinglet ~g 2 are the 



following (I follow the notation and conventions of Refs.[l7|[18| and omit the flavour 
matrices A,): 



in -i 



n L 

n-1 

- ^lf,l m D {a D^ ■ ■ ■ pw-ipw+i . . . D^-^ip] - (traces), (3) 



i=i 



= ^(Vl-Vn-l-l + ty + Un-!-!), (I = 1, ■ • • , 71 ~ 2) (4) 

= i^mS'i,^ ■ . . D^-iy^ip - (traces), (5) 

fi — 1 

ft - m)757 CT L' Ml • • • D^Y*"- 1 ^] - (traces), (6) 

where { } means complete symmetrization over the Lorentz indices and m rep- 
resents the quark mass. The symbol 5" denotes symmetrization on the indices 
[i\[i 2 ■ ■ ■ /i n -i and antisymmetrization on cr/ij. The operators in Eq.(f|) contain the 
gluon field strength and its dual tensor = \e^ va pG al3 and they are given by 

Vi = -i n gS'^ 5 D^ ■ ■■G™ ■ .•£>"»- V" -1 ^ - (traces), (7) 
U t = i n - l gS'^D^ ■ ■ ■ • • • D^- 2 ^- 1 ^ - (traces), (8) 
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where g is the QCD coupling constant. The operator R E in Eq.(^) is proportional 
to the equation of motion (EOM operator). The above twist-3 operators are not all 
independent but they are constrained by the following relation [J/J |12| : 



n-1 "- 2 



R*/ 1 -^- 1 = iC 1 " ^ 1 + E( n " 1 " O^r 1 ""^ 1 + RT'"^' 1 - (9) 



n 1=1 



Thus in total there are n independent operators contributing to the (n-l)-th moment 
of ~g 2 - But we will see later that in the Nq — > oo limit the (n-l)-th moment is 
expressed in terms of one operator R F Ul '" >ln ~ 1 . 

In all the analyses of ~g 2 performed so far in the framework of OPE and RG, 
operators Ri,R m ,RE of Eqs.@-@ have been taken as independent bases. In this 
paper I choose Rf,Ri,Re as independent operators, replacing R m with R F of Eq.([3|). 
The advantage of this choice of operator basis is that the coefficient functions take 



simple forms at the tree-level. In fact we have 17] 



££(tree) = 1, Entree) = for I = 1, • • • , n - 2 (10) 

since the anti-symmetric part of the short distance expansion for the product of two 
electromagnetic currents can be written at the tree level as 

i I d XC ^ X'(t/^ t (ic) J v (0) ) | anti— symmetric 



(11) 



where dots • ■ ■ stands for non-leading terms and 

R wi-» n -i = z -"-i^ T57 ^£)Mi . . . _ (traces) (12) 

are twist-2 operators which contribute to g\ and g%"' 2 . It is true that due to the 
relation, Eq.(|S|), i?^ w *" Mn ~ 1 can be expressed in terms of other operators. When 
eliminating R F , we obtain a different set of coefficient functions. In other words, the 
(tree-level) coefficient functions are dependent upon the choice of the independent 
operators flTfl. 
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The renormalization constants for this new set of independent operators are 
written in the matrix form as 







( ZpF 




Zfe \ 


(R n F 


R? 


- 


ZlF 


Zij 




R] 


We) 




\ o 





Zee J 


\ R% 



l,...,n-2 



(13) 



El R 



where the suffix R(B) denotes renormalized (bare) quantities. 



Now we proceed to the moment sum rule for g 2 . Define the matrix elements of 
composite operators between nucleon states with momentum p and spin s by 

(14) 
(15) 
(16) 

Normalization is such that for free quark target we have d n = 1 and /„ = 0(g 2 ). It 
is recalled that physical matrix elements of the EOM operators vanish |[20|| . Using 
Eqs.flU]) - (|16|), we can write down the moment sum rule for g 2 as, 



(p,s 


K F 




71 — 1 

d n (s a p^ s^p^p^-'-p^- 1 

n 


(p,s 


R a>J ' 1 


'^M = 


~f l (s cr p f11 — s^ il p CF )p tJ ' 2 ■ ■ ■ p^ n ^ x 


(p,s 






0. 



M n = / dxx n ~ x g 2 {x,Q 



n 



n-2 



^d n EUQ 2 ) + lj2f l n EnQ 2 ). (17) 



The coefficient functions E F (Q 2 ) and E™(Q 2 ) satisfy the following renormaliza- 
tion group equation, 

d 



^ +m Wg- lmi ~ g)m d^ 



Ei = 'jjiEj for i,j = F, 1, 



n 



where (3(g) and 7 m (5 l ) are the QCD j3 function and the anomalous dimension of 
mass operator, respectively. The anomalous dimension matrix 7^ of the composite 
operators R F and Rf with I = 1, ■ ■ ■ ,n — 2 is defined as 



~ x dz 



for i,j — F,l,---,n — 2 . 



(19) 



Note that the anomalous dimension matrix which appears in Eq.(|T8D is a transposed 
one. This comes from our convention of defining renormalization constants and 
anomalous dimensions of the operators in Eqs.(|TJ|) and flU 
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In the leading-logarithmic approximation, the solutions of the RG equations in 
Eq. ([l8|) are given as follows [pIR : 



exp 



7 



(0)n 



2/?o 



-In 



'a(Q 2 
a(fi 2 



for i 



FA, 



,n-2. 



(20) 



Fi 



where a(Q 2 ) is the QCD running coupling constant, (3q and ^°> n are, respectively, 
one-loop coefficients of the j3 function and anomalous dimension matrix, 

1 llN c -2n f 



% (g) 



O 



(47T) 



Til 



(0) V + O(g') 



(21) 
(22) 



with rif being the number of flavours, and we have used the fact that E F ({i 2 ) = 1 
and Ef(ii 2 ) = (for I — 1, • • • , n — 2) at the lowest-order. 

Now we need the information on the anomalous dimensions (^°' n )Fi (for i = 
F, 1, • • • , n—2). We can get it without embarking on a new calculation of the relevant 
Feynman diagrams. We utilize the existing results on the anomalous dimension 
matrix for the operators Ri,R m and Re- In the case of the conventional choice of 
Rk,R m and Re as independent operators, the renormalization constant matrix takes 
a triangular form 



( R l\ 


(z lj 






( R l 







7 


° 


pn 

- n 'm 


\ R e) B 


V o 





Zee ) 


\R n E 



(I, J 



n-2) 



(23) 



R 



In the MS renormalization scheme is expressed as 

,2 



3 16vr% j 



, n — 2, m, E) 



(24) 



where e — (4 — d)/2 with d the space-time dimension, and the components Xij have 
been calculated || [0 |L8| . The following is the result on taken from Ref . |L8 



X, 



G 



(j + W + 2) 
(Z + l)(/ + 2)(/-j) 

l+j n-iC-^x (n-l + l-j) 



+ (2c F - a 



G) 



-r 



+ 



2(-iy 



(n-l)(Z-j) /(/ + !)(/ + 2) 



iCj 



5 



Xn 



a 



G 



(1 < J < / - 1), 
11 1 1 Q q \ 



(25) 



+Cf (3 — 2Si — 2S n -i-i) , 
(n-l-j)(n-j) 



X,. 



(26) 



G* 



(n-l-0(n-0O'-0 



+ (2C F -C G ) 



4C, 



:-i) 



« +j n-2Cj jn-l-l+j) 



-l) n ~ 



n-2-lCn-2-j 



n-2.Ci (n - - I) n-l 
{I + 1 < j < n - 2) 



irn 



X„ 



-±c F s n . 



(27) 
(28) 



nl(l + l)(l + 2) ' 

If we impose that the renormalized and bare operators respectively satisfy the 
constraint Eq.@, we find from Eqs.fll3|) and (|23|) that Z's are related to the con- 
ventional Z's as follows: 

n n ~ 2 

Zff = Z mm -\ > . (n — 1 — l)Zi m , 

n-l f-f 



n-2 



-(rt - 1 -j)Z FF + J2( n ~ 1 ~ l ) z iji 
n „ 



i=i 



n — 1 



n 



(n-l-j)Z, 



Imi 



(29) 
(30) 

(31) 
(32) 



n — 1 

where Z, j = 1, • ■ ■ , n — 2. Using the MS scheme-rule, l/e — > ln/i 2 , we obtain from 
Eqs.fllD and (||), 

n n - 2 

X mm ^ V(n - 1 - 0-Xim> (33) 

n-l ^ 

n-2 



^ 2 7ff 



z=i 



o 2 ~(0)n 

-8tt 



5 2 ~(0)n 

^ Tip 



-in -l-j) 



X„ 



Tl 

+ — -X>-l-Z)X, m 



n-2 



+ £(n-l-Z)X y , 



i=i 
n 

n — 1 



-X, 



Irai 



(34) 
(35) 
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3 2 ~(0)n 



n 



Xij An - 1 - j)X lm , 

n — 1 



(36) 



It is straightforward to calculate the above using the expressions JQ,- in Eqs.( 



(|28 ). Especially we obtain 



^ 2 7^ 



3 2 ~(0)n 
]7T 'Fj 



(37) 



■(2Cj, - Co) 



(n-1- j)\ 2S' n _i - Sj - S^i-i-i + 1 + 



j'-i 



n-1 j(j + 1 )(j + 2) n-j 



n-2 



+ £ (n - l - ojenw-^-'^-t + '-j) 



+ ... 2( : 1)J .. , c 



z=j+i 



. 2 Q_! (n-l)(l-j) 1(1 + + 2) 



for j = 1, . . . ,n — 2 



(38) 



Now we see that the mixing anomalous dimension 7pj n turns out to be proporti- 
nal to (2Cp — Cq). Since 



2iV, 



C G = N c , 



(39) 



we have 2Cp = Co and thus 7^]™ = in the iV^ — ► 00 limit. Then Eq. (pO|) gives, 



E n F (Q 2 ) 
E?(Q 2 ) 



am 

a(/i 2 ) 







for I = 1, • • • , n — 2 . 



(40) 
(41) 



Returning to Eq.flPTD, we find that, at Nc going to infinity, the moment sum rule 
for g 2 takes a simple form as follows: 



/ dxx n 1 g 2 {x-,Q 
Jo 



' 2n r ' 



am 

a{fi 2 ) 



(42) 
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with 



~(0)n 

Iff 


2iV c 


Sn-1 


1 + 


i " 

2n 


2(3 


1 

3 


UN C - 


2n/ 





(43) 



In other words, at large Nc, the operators fij^ 1 '"^- 1 involving the gluon fields 
decouple from the evolution equation of g 2 and the whole contribution is represented 
by one type of operators Rp 11 Mn_1 . With the substitution Cf = Nq/2 and n = j + 
the anomalous dimension 87r 2 7^ n coincides with Eq.(18) of Ref.[|]j|. This completes 
the reproduction, in the framework of OPE and RG, of the ABH result on g 2 . 

It should be emphasized that we have reproduced the ABH result without as- 
suming massless quarks. A question expected to come up immediately is that the 
replacement of the mass- dependent operator R^ with R F may be equivalent to 
working with massless quarks. The answer is no. Indeed it can be shown that even 
when we include the mass-dependent operator R^ among the independent operator 
bases we reach the same conclusion. Let us take, for an example, R F , Rf (with 
I = 2, • • • , n — 2), R 7 ^ and R% as independent operators replacing one quark-gluon 
operator Rf = i with R F . With this choice of new operator bases, the moment sum 
rule for ~g 2 is written in terms of the coefficient functions E F (Q 2 ), E™(Q 2 ) with 
/ = 2, • • • , n — 2, and E^(Q 2 ). The renormalization constants for these operators 
are written as 



(R F \ 




( 




Zj? m 


Zfe \ 


(R n F 


\ 




R? 




ZlF 


Zij 


Zlm 


Z\ E 


Tin 

j 
















Zmm 


JD 


Tin 
- n 'm 






\R n E ) 


B 


I o 








Zee / 


\R n E 


J 





(l,j = 2,--.,n-2). (44) 



Again imposing that the renormalized and bare operators respectively satisfy the 
constraint Eq.(|9|), we find that Z's are related to conventional Z's as follows: 

-1 n-2 



jpF 



n-2 



5>-i-03 



i=i 



n-2 



(45) 



Z Fj = -(n-i-j)Z FF + y £(n-l-l)Zi3> (j = 2,---,n-2) (46) 

i=i 

n — 1 - n - 1 n_2 
2f m = Zpp H Z mm + } (n — 1 — l)Zi m (47) 

n n ti 
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Then it is easy to obtain the following one-loop coefficients of the relevant anomalous 
dimensions 



j 2 ~(0)n 
to 1f > 

2 2 ~(0)n 



*c F (s n . 



1 1 



4 2n> 

+ terms proportional to {2C F — Cg) 



oc (2C F -C G ) 
oc (2C F -C G ) 



for j 



,n-2 



(48) 
(49) 
(50) 



Inserting these anomalous dimensions to the solutions of the RG equations for the 
coefficient functions E F (Q 2 ), E?(Q 2 ) {I = 2, • ■ • , n - 2) and E^(Q 2 ) , 



exp 



7 



(0)n 



2/3o 



-In 



a(Q 2 



Ct(/i 2 



for i = E 2, 



2, m 



we obtain in the large N G limit 
E n F {Q 2 ) = 

e?(Q 2 ) = 
K(Q 2 ) = 



«(Q 2 




0. 



for / 



E F (Q 2 ) 
2,- 



n 



(51) 

(52) 

(53) 
(54) 



Thus we reach the same conclusion Eq. (|42|) even when we include the mass-dependent 
operators among the independent operator bases. 

A few comments are in order. Firstly, the twist-3 quark-gluon operators 
decouple from the evolution equation for ~g 2 a t large Nq. This might be explained 
by an argument on quark condensate 0. A hint is that the mixing anomalous 
dimensions 7^j n turn out to be proportinal to (2C F — C G )- There are two types in the 
products of colour matrices entering into the calculation of anomalous dimensions 
for the flavour nonsinglet ~g 2 '- 



rpbrparpb 



(C F - \c G )T a 



1 



C F T a 



^4 



2N C 

rpa 



(55) 
(56) 



2V N c < 

It is argued in Ref. [0] that the quark condensate contains all colours and at large 
N G the condensate polarization becomes small and that the combination T b T a T b is 
connected with condensate polarization effects. 
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Secondly, we have chozen particular sets of the independent operators and reached 
a simple form for the moments of g 2 in the large N c limit. However, arbitrariness 



in the choice of the operator bases should not enter into physical quantities [17 



A different choice of the operator bases leads to different forms for the anomalous 
dimension matrix and the coefficient functions. Recall that the constraint, Eq.(^), 
gives a relation among the tree-level coefficient functions and also a relation among 
the matrix elements of the operators. After diagonalizing the anomalous dimen- 
sion matrix and using these relations, we can arrive at the same conclusion for the 
moments of ~g 2 111 the Nc — > oo limit. What we did in this paper is that we choze 
particular sets of bases from the beginning which include an operator that represents 
the whole contribution to g 2 for large Nc. 

Finally, the nucleon has other twist-3 distributions, namely, chiral-odd distribu- 
tions Q 2 ) and e(x, Q 2 ) |22[. Just like the g 2 case, the Q 2 evolutions of flavour 
nonsinglet h^x^Q 2 ) and e(x,Q 2 ) turn out to be quite complicated due to mixing 
with quark-gluon operators, the number of which increases with spin. However, it 



has been proved |23j that in the large Nc limit these twist-3 distributions also obey 
a simple DGLAP equation. The proof holds true only when we work with massless 
quarks. 



This work was inspired by an interesting talk given by Y. Koike at International 
Symposium on QCD Corrections and New Physics, 27-29 October 1997, Hiroshima. 
I would like to thank him and also thank the organizer of the Symposium, J. Kodaira. 
The discussions with Y. Koike and T. Uematsu on the twist-3 operators in the 
large Nc limit were indispensable for the completion of this paper and are happily 
acknowledged. This work is supported in part by the Monbusho Grant-in-Aid for 
Scientific Research No. (C)(2)-09640342. 
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